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Mullins has treated a two-dimensional model of grain boundary motion in which each point on the 
boundary moves toward its center of curvature with a speed proportional to its curvature. For bound- 
aries which preserve shape under uniform magnification, an integral representation of the boundary 
shape is found. We then obtain several analytic results from approximate evaluations of the integral. 
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The motion of an idealized grain boundary in two 
dimensions has been treated by Mullins fl]. 1 In this 
model, a given point on the grain boundary moves 
towards its center of curvature with a speed that is 
proportional to the curvature. Recently Sun and Bauer 
[2] have extended the Mullins treatment by numeri- 
cally solving the differential equation for the boundary 
shape for the case of boundaries which are invariant 
under magnification. They have used the results to 
measure grain boundary mobilities [3], 

The shape of a boundary which is invariant under 
magnification is given by 



r(0, t) = R(0) T(t) 



(1) 



where r and 6 are polar coordinates and t is the time. 
T(t) = |T 2 (0) + 2ct] 11 ' 2 , where c is a separation constant 
and R(6) satisfies the nonlinear differential equation 



{R* + 2RI-RRee}/{R 2 [R 2 + R 2 e ]}=-(clk). 



(2) 



Here, the subscript 6 denotes differentiation with 
respect to 6 and k = Mcr where M is a mobility and o~ 
is the boundary free energy. 

The purpose of this note is to obtain some analytic 
results for R(0); in particular, the solution of eq (2) 
can be reduced to a single quadrature. Following 
standard techniques [4], we find 



Jr 



0(R)-6o= da a~ l {C,a 



exp(ca 2 /k) — 1} 



1/2 



(3) 



where 0(/?o) = #<> and C\ is an integration constant. 
Sun and Bauer consider that a free surface acts as a 
mirror plane; thus for a free surface along the ray 
= do we take as an additional boundary condition 
Re = at Ro. Introducing the change of variable 

7] = (a' 2 /RD — 1, eq (3) becomes 

0(R)-0o=i [dT7(l + T7)- ! {(l+*?) 
J o 



exp (br)) — 1) 



-1/2 



(4) 



1 Figures in brackets indicate the literature references at the end of this paper. 



where p= {H' 2 IR 2 ) - 1 and b = (elk) R 2 . The boundary 
shapes for different values of b have been discussed by 
Mullins. Figure 3 of reference [1] and figure 1 of refer- 
ence [2] illustrate boundary shapes for various values of 
b. In reference [1], Ro = S/2 and 00 = 77-/2 and in refer- 
ence [2], 0o = O and a and f(a) are equal to our 0(oo) — O 
and b, respectively. From the equation for T(t), we see 
that b > yields solutions which increase in size with 
time while b < gives solutions which shrink toward 
the origin r=0. For example, from eq. (2) it is clear that 
6 =— 1 is a circle of radius /? ; eq (4) is singular in this 
case since when R(6) =R , * 0(R). The case b = 
corresponds to a straight line. For b > 0, Sun and 
Bauer [2] solved numerically a differential equation, 
equivalent to eq (2), to obtain the boundary shape. For 
b > 0, R ranges from R {) to infinity. For b < 0, R ranges 
from Ro to a value of R at which the quantity [(1 +p) 
exp(6p) — 1] becomes negative. It follows that for 
— 1 < b < 0, R is greater than R while for b < — 1, R is 
less than R G . 

Although it does not appear possible to evaluate the 
right-hand side of eq (4) analytically, we can obtain 
some approximate analytic results. For values \bp\ 
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=\(blRh)(R- 2 -Ri)\«l(and 6^-1), we have (1 + tj) 
exp (fr»j)-ls (1 + 6)7? {H-(6tj) (l + 6/2)/(l + 6)} 
and integration of eq (4) gives 



cos {2(l-6 2 /2)" 2 (6>-6>o)} 

= {l-p(l-6-6 2 )/(l + 6)}/{l+p}- 



(5) 



When 6 » 1, the only contribution to the integrand 
occurs for r) near zero. We can than approximate the 
integrand by {exp (br\) — 1} _1/2 . Upon integration and 
rearrangement, we obtain 

R=Ro{l-(2/b) In cos [6(0- 0o)]} 1/2 (6) 

We see that 6 — 6o=7rl2b when R is infinite. Thus, 
eq (6) gives the shape of a grain boundary which makes 
a very small angle with a free surface. 

Finally, for R infinite, we obtain an approximate 
analytic formula for 0(°°) — do as a function of b. 
We rewrite the integrand of eq (4) in the form 
[(1 + b)^]-'! 2 [1 + t?]- 1 exp [- bj)l2] [1 + {1 - bi) 
— e- b7) }l{i){\ + 6)}] _1/2 . Expanding the last factor in 
powers of br\, yields 

[l + {l-67/-e-^}/{7?(l + 6)}]- 1/2 

= 1 + q:iT7 + q:2T7 2 + . . ., 

where a! = 6 2 /4(l + 6) and a 2 = {6 3 /96(l + 6)}{[96/ 
(1 + 6)] —8}. The resulting integrals can be evaluated 
in terms of error functions to give 

eM-6o={7rl2b(l + b)}*IHg+a 1 (l-g) 

+ a 2 te-1 + 6- 1 )}, (7) 



where g= (tt-6/2) 1 / 2 exp (6/2) erfc [(6/2) 1 / 2 ]. The a 2 
term is small compared to the leading term. Thus, eq (7) 
provides an analytic representation of the numerical 
results for b as a function of 0(°°) _ #o, which are given 
in figure 2 of reference [2]. It is perhaps surprising that 
eq (7) is a good approximation for all values of 6, since 
it is based on an expansion in powers of (brj). However, 
when brj is not small, the integrand is very small due to 
the exponential term and makes a negligible contribu- 
tion to the integral. For b < < 1 , eq (7) yields 



0(cc)-e o =j [1-(26/tt) 1 / 2 ] 

and for 

6^oo, 0(oo)-0 o =&-i 

(77/2) 1 / 2 [l+ (1/4) + (1/32)]= 1.61/6, 



which is in reasonable agreement with the result 
(7r/26) obtained from eq (6). As an example of an inter- 
mediate value of 6, we find that for 6 = 7.1, 0(o°) — O = 
10.0° in agreement with Sun and Bauer. 
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